Generalized Potentials for a Mean-field Density Functional Theory of a Three-Phase 

Contact Line 
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We investigate generalized potentials for a mean-field density functional theory of a three-phase 
contact line. Compared to the symmetrical potential introduced in our previous article p[, the three 
minima of these potentials form a small triangle located arbitrarily within the Gibbs triangle, which 
is more realistic for ternary fluid systems. We multiply linear functions that vanish at edges and 
vertices of the small triangle, yielding potentials in the form of quartic polynomials. We find that 
a subset of such potentials has simple analytic far-field solutions, and is a linear transformation of 
our original potential. By scaling, we can relate their solutions to those of our original potential. 
For special cases, the lengths of the sides of the small triangle are proportional to the corresponding 
interfacial tensions. For the case of equal interfacial tensions, we calculate a line tension that is 
proportional to the area of the small triangle. 
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I. INTRODUCTION 

In our previous article [TJ, we treated a three-phase 
contact line by a mean-field density functional that in- 
volves a symmetric potential. This potential allows us 
to find analytical far-field solutions because one of the 
mole fractions is constant for any two-phase transition. 
In order to develop a phenomenological model for a non- 
symmetric potential, which is more realistic for ternary 
fluid systems, we need a more general potential that has 
three minima in arbitrary locations within the Gibbs tri- 
angle. In this article, we introduce a flexible approach to 
construct such generalized potentials. 

As illustrated geometrically in Fig. 



1(a) a three-phase 



contact line is modeled for a ternary fluid system having 
three bulk phases a, /3, and 7, which subtend dihedral 
angles 6 a , Op, and Q 1 as sketched in Fig. l(a)| Each of 
the interfaces a/3, f3j and 7a, if extrapolated far from 
the contact line, is perpendicular to a side of the dashed 
triangle, which is illustrative of our much larger actual 
computational domain. The contact line is where the 
three interfaces appear to meet. We assume the system 
to be translationally invariant along the vertical direction 
perpendicular to this figure, thus reducing this problem 
to two dimensions. Compared to the homogeneous bulk 
phases, the inhomogeneity arises from the formation of 
interfaces and a contact line, which are actually diffuse 
regions. 

Based on the thermodynamic method introduced by 
Gibbs @, p. 228], the inhomogeneity is treated by means 
of an excess grand potential £l X s> which, by convention 
@, ch. 8], can be expressed by 

Qxs = Lt + LR a pa a /3 + LRp 1 ap 1 + LR ja a ja . (1) 

Here, L is the length of the contact line and Rij is the 
distance from the contact line along interface ij toward 



the boundary. In the limit of all Ri 



00, the line 



tension r is defined as the excess grand potential per 
unit length associated with the contact line while each of 
the three interfacial tensions (excess grand potentials 
per unit area) is associated with an interface ij in the 
far- field limit. A classical result [1] shows that 
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According to this result, the boundary is actually a Neu- 
mann triangle since the three sides are proportional to 
three interfacial tensions. 

Besides extensive studies of interfacial tensions (See 
Q), the line tension of a contact line plays a crucial role 
in a broad range of physical phenomena such as the equi- 
librium shapes of small droplets [1,0] > microfluidics [1,11], 
heterogeneous nucleation [l(|, cell adhesion [ll|, the dy- 
namics of a drop spreading on a liquid thin film [I2T Il3j 
and the behavior of line tension at wetting transitions 
[T3 |. For a review for both experimental and theoreti- 
cal aspects of line tension, see [l5| , and for a conceptual 
review, see 1161. 



A. Mean-field Density Functional Theory 

According to the mean-field density functional theory 
introduced in our earlier work [l|, the excess grand po- 
tential for a ternary fluid system is 



a, 



BL [/ + g] dA, 



(3) 



where / is a potential function and g is a gradient energy 
of the chemical constituents. In this article, we assumed 
that the dominant intermolecular forces in our system 
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FIG. 1. Schematic diagram for the mapping between phys- 
ical space and the Gibbs triangle (briefly, the Gibbs space), 
(a) Sketch of our system consisting of three bulk phases a, f3, 
and 7, divided by three interfaces a/3, /?7 and 7a. The inter- 
faces, if extrapolated, are each perpendicular to sides of the 
dashed triangle (illustrative computational domain) and meet 
at a three-phase contact line. The corresponding dihedral an- 
gles are 8 a , 6/3, and 9~, and (x,y) are Cartesian coordinates 
in physical space. The region within three solid curves rep- 
resents approxi matel y the diffuse region of the interfaces and 
the contact line, (b) The Gibbs space, in which the three bulk 
phases are located at three points and (u, v) are Cartesian co- 
ordinates in units of mole fractions. The three curved lines 
that connect pairs of points for bulk phases are trajectories of 
two-phase transitions, (c) Slice of an interface in a cuboid far 



from the contact line, s is a coordinate perpendicular to the 
a/3-interface and w is the width of the cuboid in a direction 
perpendicular to the contact line. 



are short range, so physical quantities can be formulated 
in terms of local densities. Based on the assumption 
of uniform molar volume, we change the variables of / 
and g from number densities to mole fractions JQ, and 
2i=i* — 1 ( see 0i sec - II] f° r details). For reviews 
of the general mean-field density functional method of 
interfaces, see |3. ITm]- 

In this phenomenological theory, there is no definite 
form of the potential / for the excess grand potential Q xs 
in ([3]). The specific form that we used in earlier work [l| 



f = J2(X i -a) 2 (X i -b) 2 = f°, 



(4) 



where < a < 1 is a constant and b = (1 — »)/2, is 
a symmetric quartic function that allowed us to obtain 
asymptotic analytical solutions in the far-field. This form 
is an extension of the potential used in the Landau's phe- 
nomenological theory [TH [l9| . Similar two-density quar- 
tic potentials have been used in [13, HH for line tension 
and for a first order wetting transition. Sixth order poly- 



nomials have been used for second and higher order tran- 
sitions [2ll - [2"3l | . In this article, we focus on potentials of 
quartic form. 

Compared to other two-density models, our model is 
actually pseudo-binary since it formulates the potential 
in terms of three mole fractions {-Xj}i=i,2,3 by the as- 
sumption of uniform molar volume. Because of the con- 
straint Yli=i — 1, there are two independent mole 
fractions, so our potential can be described in terms of 
two independent variables, although we sometimes dis- 
play all three mole factions to illustrate its symmetry. 

Note that the gradient energy g in the functional (|3]) 
of the excess grand potential Q xs is specified by 



(5) 



where {^j}j=i,2,3 are constants. For cases with equal 1% 
(isotropic gradient energy) , the three- fold symmetric po- 
tential leads to three-fold symmetry of the physical do- 
main. To resolve this special geometry, we employed a 
triangular grid to obtain numerical solutions over the en- 
tire domain (for details, see [H, sec. III]). 



B. Mapping from Physical Space to the Gibbs 
Space 

In our mean-field density functional theory of a three- 
phase contact line for a ternary fluid system (summarized 
in Sec. II A( for details, see [H, sec. II] ), the functional of 
the excess grand potential consists of a potential function 
and a gradient energy. Since we use three mole fractions 
as variables, there is a connection between the physical 
space and the space of the Gibbs triangle (briefly, the 
Gibbs space). This connection depends on the specific 
choice of potential and the coefficients of the gradient 
energy. Fig. Q] illustrates a mapping between the physi- 



cal space of our system 1(a) and the Gibbs space 1(b) 



In general, the three bulk phases of the physical system 
are represented by the three minima of a given poten- 
tial, which are located at three points within the Gibbs 
triangle. 



The curved lines within the Gibbs triangle (Fig. 1(b) ) 



are the trajectories of the two phase-transitions in the 
far-field regime between any pair of bulk phases in the 
physical space (Fig. 1(a) ). The far-field regime is located 
at a distance that is far from the contact line compared 
to the interfacial widths (or the size of the central core 
region associated to the contact line). In the far-field 
regime, a slice of interface can be contained in a cuboid 
(Fig. 1(c) I and the interfacial width is constant. For a 



two-phase transition between any pair of bulk phases, the 
change of mole fractions from one bulk phase to another 
maps to a curved line connecting two potential minima 
within the Gibbs triangle. At equilibrium, these trajecto- 
ries minimize the excess grand potential according to the 
form of the potential and the coefficients of the gradient 
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energy. Since the interfacial tensions are excess grand 
potentials per unit area, their values depend on these 
trajectories. 

According to the form ([1]) of the excess grand potential, 
the line tension is the residue of excess potential in which 
we subtract the contribution from the far-field interfacial 



tensions. In physical space (Fig. 1(a) I, it relates to the 



core of the diffuse region centered at the contact line 
as if we subtract the cuboids of the interfaces extended 
from the far-field until they meet at the "contact line". 
This core region corresponds to three-phase transitions 
among all three bulk phases in the physical space. It 
can be mapped to a region within the Gibbs triangle 
(Fig. 1(b)) surrounded by the three trajectories of the 
two-phase transitions. Similar to the interfacial tensions, 
the line tension is the excess grand potential per unit 
length. Its value depends on the form of the potential 
within this core area and, of course, the coefficients of 
gradient energy. 

In this article, we extend our previous model to more 
realistic systems. Specifically, we use geometrical reason- 
ing to construct generalized quartic potentials with three 
minima arbitrarily located within the Gibbs triangle. For 
a subset of these potentials, the resulting potential is a 
linear transformation of our original potential and we 
can obtain simple analytic far-field solutions. We con- 
nect these solutions to our original potential by scaling. 
For some special cases, we relate interfacial tensions and 
line tension to the lengths of the sides and the area of 
the small triangle formed by the three minima. 



II. GENERALIZED QUARTIC POTENTIALS 

To find more general quartic potentials with arbitrary 
minima, while maintaining the condition of constant mo- 
lar density, we start with a discussion of quartic functions 
with two independent variables that correspond to two 
independent mole fractions. In general, a two-variable 
quartic polynomial can be expressed by 



f(x,Y) = j2J2 a ^-* xtYn ~ 



(0) 



n=0 i=0 



which contains 15 independent parameters ai^ n -i. How- 
ever, we need a form such that / is positive except 
for / = at the three minima, which gives us three 
equations relating the coefficients. At each of the 
minima, we require df/dX — df/dY — 0, which 
amounts to six conditions. Furthermore, in order to 
have parabolic potential wells, we need the second deriva- 
tives of / to satisfy the following inequalities at the 
three minima: d 2 f/dX 2 > (or d 2 f/dY 2 > 0) and 
(d 2 f/dX 2 )(d 2 f/dY 2 ) - (d 2 f/dXdY) 2 > 0. This gives 
us six inequalities. In general, we have 15 — 3 — 6 = 6 
free parameters together with six additional constraints. 

In this framework, our original potential in Q is a 
special case which has only one parameter, a, represent- 
ing the size or orientation (magnification or inversion) of 



the equilateral triangle formed by the three minima 0, 
sec. II. A, discussion after eq. 14]. In the development to 
follow, we obtain a positive potential by assuming it to 
be a sum of squares of various expressions. We break the 
symmetry of our original potential by locating the three 
minima at the vertices of a small triangle having any 
shape and orientation within the Gibbs triangle. In the 
following discussion, we use the term a "small triangle" 
to denote an inner triangle formed by the three minima 
of a given potential within the Gibbs triangle. 



A. First Generalization 

At first, we explore the structure of our original poten- 
tial f° in (Q| . Xi can be expressed as a function of two in- 
dependent Cartesian variables u and v, i.e. Xi — Xi(u, v) 
as illustrated in Fig. 1(b) In [l|, we showed that f° can 
be scaled to the form 



/ = 



f" 



(a -by 



J2Y i (u,v) 2 (Y(u,v)-l) 



(7) 



where Yi(u,v) — (Xi(u,v) — b)/{a — b) are scaled mole 
fractions, and X^ =1 li = 1. This form is equivalent to 
f° when a — 1. For simplicity, we mainly compare new 
potentials to /. 



As illustrated in Fig. 2(a) we specify explicitly the 
Cartesian coordinates (u, v) for 5^ (it, v) 



y/3 v 
Yi{u,v) = — —u - - + 1 

Y 2 {u,v) = —u - g 
Y 3 (u,v) = v. 



(8) 



In l[8]l, the symbols u and v bear the same relationship 
to the scaled mole fractions Yi as u and v in Fig 1(b) 
do to the unsealed mole fractions X t . The three minima 
of the potential / in Q, namely a = (u^vf) = (0,0), 
13 = (u° 2 ,v° 2 ) = (2/^3,0), and 7 = (u%v%) = (1/^3,1) 
are located at the three corners of the Gibbs triangle as 
in Fig. 2(a) consistent with the contour plot in Fig. |2(b)] 
Note that *i("°> w °) = Y i(uli v k) = for 3 k ^ h and 
Yi(v$,i#) = l. 

If we define two sets of linear functions L®(u,v) = 
Yi(u,v) and Z?(u,u) = Y,(u,v) - 1 = L?(u,v) - 1, the 
form ([7]) of the potential / can be expressed as 



;=$>°M) 2 /°( 



u,v) 



(9) 



As shown in Fig. 2(a) L\ — is a line in the Gibbs space 
that passes through two minima (upV®) and (it^,u°), 
where j ' ^ k ^ i, and if — is a line parallel to Lj = 
that passes through the remaining minimum (it?, v®). For 
each minimum of /, there are three intersecting lines. 
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(a) Geometry of the three minima of / (or f° when 
a = 1) 




(b) Contours of / 

FIG. 2. | (a) | Geometry of the three minima of the potential / 
given by @ (or f° in Q when a = 1). / is a function of scaled 
mole fractions {K;};=i,2,3, where Y^ =1 Yi = 1. Yi can be 
expressed in terms of two independent Cartesian coordinates 
(«,«). a = («?,«?) = (0,0), 13 = {u° 2 ,v° 2 ) = (2/V3.0), and 
7 = (w3,«3) = (l/-\/3, 1) are the three minima of / located 
at the three corners of the Gibbs triangle. Li — represents 
a line that passes through two points (m°,v°) and (u°,v k ) 
(j k i), and L — is a line parallel to = but passing 
through only one point (Ui,v^). (b) Contour plot of /. The 



value of / decreases from its value at the center of the inner 
area as one proceeds toward the three vertices where it is zero. 



Based on the geometrical interpretation of the form © 
of the potential /, we can generalize its structure by al- 
lowing these pairs of parallel lines to move and requiring 
the intersections of three lines to be located within the 
Gibbs triangle. To generate a potential function with 
the desired properties, we locate the three minima ar- 
bitrarily, a = (ui,vi), (3 = (u 2 ,v 2 ), and 7 = (1*3,^3), 



within the Gibbs triangle the potential . Then, we define 
two sets of linear functions Li(u,v) and Ii(u,v), where 
Li{u, v) = is a line passing through two minima (v,j, Vj) 
and (life, Ufc) for j ^ k ^ i, and Ii(u, v) — is a line par- 
allel to Li{u,v) — and passing through the remaining 
minimum (ui, Vi). 

As illustrated in Fig. 3(a)| the line Li — coincides 
with the side of the small triangle opposite to the vertex 
i, while the line Li — passes through the vertex i and is 
parallel to the line Li = 0. Thus, a generalized potential 
with arbitrary mimina is given by 



fg = ^2 dfLi(u, v) 2 Ii(u, v) 2 



(10) 



»=i 



where the di are nonzero weighting coefficients that re- 
late to the curvatures along the lines that connect each 
pair of potential wells. For each term, Li(u, v) Ii[u, v) 
represents a positive quartic function that vanishes along 
two parallel lines Li(u, v ) = and Ii(u, v) = 0. We need 
all three terms to produce a positive quartic function that 
vanishes at the three points (it;, 1*1)1=1,2,3 which are the 
intersections of three lines chosen from the three pairs 
of parallel lines. Contours for such a potential are illus- 



trated in Fig. 3(b) 



The choice of Ll(u,v) and Ii(u,v) for the potential f g 
in (|10[) is somewhat arbitrary but we use specific forms 
of Li(u,v) and Li(u,v) and a weighting coefficient for 
each term to maintain generality. The specific forms of 
Li(u,v) and Ii(u,v) that we use are 



Li(u,v) ee Cijk[{u - Uj)vjk — {v- Vj)u 3k } 
= C ijk [(u - u k )vj k - (v - v k )uj k ] 

I%(u,v) ee Li(u,v) - h 

= C ijk [(u - Ui)v jk - (v - Vi)ujk], 



(11) 



where C ijk = (e 

ijk "h I £%jk 

|)/2, eijk is the Levi-Civita sym- 
bol for three dimensions, Uj k ee iij — u k and Vj k = Vj — v k ; 
moreover, h = Lj(u,-,u,-) = 2A, and A is the area of the 
small triangle. For simplicity, we require the indices of 
the three vertices to be numbered counterclockwise. Note 
that 



h = ^2 £ijkUjV k 
(i,j,fc)e{i,2,3} 



1 1 1 

Ui u 2 u 3 
Vi v 2 v 3 



(12) 



If we express Li and h in terms of two independent 
mole fractions Y\ and Y 2l we obtain 



Li = -j=C ijk [(Yi - Yl k) )Y 2 W - (Y 2 Y 2 {k) )Yl 



VJ (i,i,fc)e{i,2,3} 



in which we have defined Y- 



(jk) 



Y 



(i) 



(13) 



X. (fe) and 



expressed the three minima (v,i, Uj)i=i 2,3 m terms of 
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(a) Geometry of the three minima of f g 




(b) Contours of f g 



(a) Geometry of the three minima of the quartic 
= (ui,vi), P 



FIG. 3. 

potential f g given by (fTO)) . a = (tti,«i), /3 = (ii2,«2), and 
7 = (v,3, V3) are the three minima of this potential arbitrarily 
located at any three points within the Gibbs triangle. Li — 
represents a line that passes through (uj,Vj) and (uk,Vk) for 



j 7^ k 7^ i, and I, = is a line parallel to L, 
ing through the vertex (u%,v, 
f g for (m,vi) 



but pass- 
(b) Contour plot of potential 
(0.3087, 0.1667j7(«2, v 2 ) 



(0.9060, 0.1167) 
(u 3 ,v 3 ) = (0.4974,0.6667), and {d u d 2 ,d 3 ) = (1.0,1.1,1.2). 



in which the three minima form an equilateral small tri- 
angle with one parameter a related to its size. With this 
parameter, we can scale the small triangle from the Gibbs 
triangle to its center point with zero area. To add degrees 
of freedom to allow this equilateral small triangle to move 
around within the Gibbs triangle, we can add two extra 
parameters corresponding to translations in two perpen- 
dicular directions. For arbitrary shape and orientation 
of the small triangle, we need a parameter correspond- 
ing to rotation and two parameters for distortion, such 
as changes of two inner angles. In general, these opera- 
tions correspond to a linear transformation with six free 
parameters, which will be discussed in Section UlI Bl 

The general properties of the potential f g in (fTUf have 
been tested by performing calculations of its first and sec- 
ond derivatives with respect to u and v. For details, see 
[24l sec. 5.1.1]. From these calculations, we prove that 
f g , which is a generalization of potential / in (J7J), sat- 
isfies our desired properties of a quartic polynomial, as 
discussed in the beginning of Sec. [TTJ We also calculated 
the eigenvalues of the Hessian matrix H, which is the 
matrix in terms of the second derivatives of f g with re- 
spect to u and v. These eigenvalues A± are the principle 
curvatures at the mimina, i.e. 



A± =h 2 


U=i 






± 









(14) 



where Si = Cijk^Juj k + v l - k is the length of the side of 

the small triangle opposite to the vertex as in 

Fig. [8] Note that both of its eigenvalues A± are positive. 

For the case that all di = 1 and S\ = S2 = S3 = 
(2/y / 3)(a — b), the eigenvalues of H in (fl~4| reduce to 



(5j ,Y% f )j=i 2 3 by using the relation |5J of (u,v) and 
Yi. 

The form (JTDJ) of the potential f g has three coefficients 
di and six parameters (it*, Vi) that locate the positions of 
the three minima, nine parameters in total. One could 
take out a common factor without changing the nature 
of the potential. From the discussion in Page [3] for the 
desired properties of a general quartic potential given by 
©, we know that there are six free parameters, which 
means there are three parameters that might be able to 
be merged with others. 

From the forms ([7]) and §§§ of the potential /, our 
original potential f° divided by a factor (a — b) A can be 
viewed as a special case of the potential f g given by (fit))) . 



A± = ^(a-6) 6 (2±V3) (15) 

which connects to our original potential f° in ([4]) . When 
a = b = 1/3, these two eigenvalues vanish, which means 
the curvatures vanish. This is equivalent to the bulk crit- 
icality for f°, for which the three minima merge to a one 
minimum located at the center of the Gibbs triangle. In 
this case, the three mole fractions are uniform through- 
out the physical space, i.e. Xi = 1/3. The resulting 
intefacial and line tensions at the vicinity of this critical 
value scale consistently with previously determined re- 
sults in our earlier work [l|, in which the ratio of their 
exponents in terms of \Xi — 1/3| satisfies the mean- field 
approximation in [25| . 
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,Z,=0 



I,=0 
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(a) Geometry of the three minima of /c 




(b) Contours of fa 



(a) Geometry of the three minima of the quartic 

= («l,Ul), /3 = (U2,V 2 ), 



FIG. 4. 

potential /g given by (fT7)> . a 
and 7 = (143,113) are the three minima of this potential ar- 
bitrarily located at any three points within the Gibbs tri 
angle. Li = represents a line passing through (uj 



(b) 



and (uk,Vk), and Ki = is a point located at (ui,v- 
Contour plot of potential fa for (ui,«i) = (0.3087,0.1667} 
(u 2 ,v 2 ) = (0.9060,0.1167), (u 3 ,v 3 ) = (0.4974,0.6667), and 
(<n,6i, Ci) = (1,1,1), (i = 1,2,3). 



B. Second Generalization 

The potential f g in (|10[) can be generalized by replacing 
Ii with an arbitrary linear function that passes through a 
point (ui,Vi), which is not necessarily parallel ot Li. For 
details, see Appendix [A] We also observed that Li(u,v) 2 
of f g does not have to vanish on a line Ii(u,v) = 0. It 
can be replaced by a positive quadratic function Ki(u, v) 
which is a paraboloid in three dimensional space that 
vanishes only at the point (iti,Vj). In general, Ki(u,v) 
can be expressed as 

Ki(u,v) = [a t {u - Ui) + b t (v - v^] 2 + c 2 (v - v{) 2 , (16) 



where a,, 6 i; and c, are three coefficients, and Cj ^ 0. 
Then a more general potential, as illustrated in Fig. 4(a)| 
can be expressed in the form 



}g = ^ v) 2 Ki(u, v). 

i=l 



(17) 



/g in (fT7|) is the most general potential we have ex- 
plored. It has nine coefficients and six parameters (it,, Vi) 
for the positions of the three minima. Thus, it has 15 co- 
efficients in total. We note that when a, = CijkdiVjk, 
hi = —CijkdiUjk, and Ci — 0, /g reduces to the potential 
f g in (JTUJ) . Therefore, / s is a subset of fc, where fa has 
six extra parameters corresponding to the orientations 
of the three potential wells, as well as the general shape 
of the potential. Recall the properties for the desired 
potential: the potential has to be positive throughout 
the entire domain and has three minima at zero poten- 
tial; at these minima, the first derivatives vanish and the 
Hessien matrices for second derivatives are positive defi- 
nite. These properties have been tested in [2J, sec. 5.1.2]. 
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FIG. 5. Schematic representation of the generalization from 
a potential / in Q to the potential f g in (|10|) of similar form 
and to the more general potential fa in (|17[) . At the LHS 
of this figure, there are three pairs of two parallel lines. The 
three minima of / are located at the three nodes where lines 
intersect and form a small equilateral triangle. In the middle 
of this figure, we generalize / to f g by allowing the three 
minima of f g to be at arbitrary locations. In this case, the 
three pairs of parallel lines determine the three minima. The 
RHS of this figure shows a further generalization in which 
each minimum of fa is located by two lines and a point. 

Fig. [5] presents a schematic summary of the general- 
ization from the quartic potential / in {7J (equivalent to 
our original potential f° in {1} when a = 1) to a simi- 
lar form (JTUJ) of the potential f g and to a more general 
quartic form p7[) of the potential fa- In the LHS of 
Fig. we see that the three minima of / form an equi- 
lateral small triangle and each minimum is located at an 
intersection of three lines. These lines consist of three 
pairs of two parallel lines. For each pair, there is a line 
passing through two minima and a parallel line passing 
through the remaining minimum. Each line corresponds 
to a paraboloid that vanishes at this line, and each pair 
of parallel lines corresponds to the minima of a quartic 
function. Thus, / is a sum of three quartic function. As 
a generalization, in the middle of Fig. [5J the three min- 
ima of f g form a small triangle of arbitrary shape. The 
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rules of selecting parallel pairs of lines and constructing 
corresponding quartic functions for f g are the same as 
for /. The RHS of Fig. [5] shows a further generalized 
potential fa, where each minimum is not located at an 
intersection of three lines but two lines and a point. In 
this case, each pair of parallel lines in f g is transformed 
to a line passing through two minima and a point at the 
remaining minimum. To compose fa, we need two sets 
of quadratic functions. The one in first set vanishes along 
a line and the one in the second set vanishes at a point. 



III. GENERALIZED POTENTIALS WITH 
ANALYTIC FAR-FIELD SOLUTIONS 

A. Straight line trajectory of two-phase transition 

In section [TTl we introduced a systematic way to find 
quartic generalized potentials by means of geometric con- 
siderations. In order to find the asymptotic solutions for 
two-phase transitions in physical space for each of the 
three interfaces in the regions far from the three-phase 
contact line, we need to know the corresponding trajec- 
tories for transitions between the pairs of wells of a given 
potential within the Gibbs triangle. According to the 
discussion in [U, sec. 2.4.1], the property of our original 
potential that leads us to simple far-field analytic solu- 
tions is that the trajectory of a transition between any 
two phases follows a straight line within the Gibbs trian- 
gle. From a geometric aspect, this straight line lies along 
a valley connecting a pair of minima; as one goes along 
this valley, the potential rises from a minimum, goes over 
a saddle point, and decreases, to another minimum. The 
first derivative and the second derivative of the poten- 
tial with respect to the normal of this straight line in the 
Gibbs space are zero and positive, respectively. For the 
subset of generalized potentials that share this common 
property, we can find analytical far-field solutions. As 
shown in Appendix [Bj the form of this subset is actually 
equivalent to the form dTU|) of the generalized potential 

u 



B. Linear transformation 

According to the geometric representation in sectionlTIl 
as illustrated in Fig. [5J the form © of the potential / 
can be generalized to a similar quartic potential f g in 
(fTUJ) with arbitrary minima («i,ui), (1x2,1*2), and (1x3, 1*3) 
by choosing proper linear functions Li(u,v) and Ii(u,v). 
In this way, the obtained potential has six parameters 
{ui, Uj}i=i,2,3 and three weighting coefficients {di}i=i,2,3, 
which match our earlier reasoning about using six pa- 
rameters that correspond to operations such as scaling, 
translation, rotation, and distortion. Furthermore, we 
notice that these operations may be represented by lin- 
ear transformations of our original potential. 



To construct a linear transformation that maps the 
three minima from the corners of the Gibbs triangle 
to three arbitrary internal points (1x1,1*1), (142,1*2)) and 
(^3,1*3) as illustrated in Fig. [SJ we write 



(18) 



where a, /J, 7, 5, £, and r\ are six undetermined pa- 
rameters. By substituting (1x1, t*i), (142,1*2), and (1*3,113) 
for (u',v'), and the corresponding three corners (0,0), 
(2^/3,0), and (1/V3.1) for (u,v), we solve 



V" 




a P 




u 






1/ 




7 5 




V 


+ 


p. 



V* 




v' 





vg , 

2 L 



\ui2 + U 31 




u 


+ 




\v\2 + V 31 




V 




yi_ 



where Uij = 
inversion, we find 

2 



Uj and Vij = Vi 



(19) 



Vj. Then, by matrix 



V3h 



5^12+^31 -fui2-«31 
Vj3„, v3 „, 
—^12 — T Ul2 



u — u\ 

V 1 — 1*1 



(20) 

where h 7^ is twice the area of the small triangle, defined 
in jnj and ([j"2]>. 



(1/V3.1) 




FIG. 6. Schematic linear transformation from the three cor- 
ners of the Gibbs triangle to three arbitrary internal points 
(ui,vi), (it2,«2), and (113,1*3) within it. 

A generalization of the potential / in j7]) with weight- 
ing coefficients a can be expressed in terms of two inde- 
pendent Cartesian mole fractions u and v as follows: 



P(u,v) = PiY^v^Y^v), y 3 («,«)) 
3 



(21) 



To construct a new potential with three minima located 
at (1x1,1*1), (1*2,^2), and (1*3,1)3), we replace Yi(u,v) in 
(HTJ by 



Zi(u',v') =Yi{u{u',v'),v{u',v')), 



(22) 



where the linear transformation (|20[) from (u' , v') to (it, v) 
has been used. Then (see 2J| for details), we find (after 
dropping the primes) 



Zi(u,v) 



Li(u,v) 



(23) 
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where, as defined in (|TT|) . Li is a linear function that 
satisfies Li(uj,Vj) = and L^u^jV^) = for (j ^ ^ 
i); whereas Lj(uj,«j) = /i, which leads to Zi(ui,Vi) = 
1. Note that Y^i=i^i = 1- Thus, the new potential is 
defined as 



f*(u,v) = f\Z 1 {u,v),Z 2 {u,v),Zj,{u,v)) 

3 

= ^ e ? Zi(U;U) 2 (Zj(W;W) _ 1) . 



(24) 



i=l 



With the definition (f2"4")l of the new potential /*, we 
find that = 0},-=i 2,3, which shows that the 

three mimina of /* are located at {(lij, i>i)}i=i,2,3- By 
expressing /*, which relates to the potential p in (l2Tj) 
(generalized from the potential / in {?])) by a linear trans- 
formation, in terms of Li, i.e. 



8=1 



^4 -n-' ^dfLili — f g , (25) 



where d; = ei/h 2 , we obtain exactly the potential f g in 
(J9j> . Thus, f g is a linear transformation from the original 
potential / with extra weighting coefficients. This is the 
reason that / and f g share a common property: the two- 
phase transitions in the physical domain follow straight 
lines in the Gibbs space. 

Notice that the structure of the two-density quartic 
potential used by Widom et al [2(J Hl[ actually belongs 
to the general category indicated by the form (|TU1) of the 
potential f g . The two relative densities of their poten- 
tial can be treated as our two Cartesian mole fractions u 
and v. We can scale their potential in terms of three di- 
mensionless quantities similar to the Zi in our expression 
and obtain a form similar to f g . However, their systems 
do not follow the assumption of uniform molar volume. 
Therefore, in their systems, the sum of the quantities 
Zi 1 and cannot be mapped into our potentials. 



C. Scaling 



From Sec. IIIIB1 we know that the potential f g in (fTU)) 
and the potential / in (J7J can be connected by a linear 
transformation. In general, a problem involved f g can 
be solved by scaling into a form like /, which we studied 
in [J sec. IV. C]. As an example, we consider the excess 
grand potential given by ([3]) but with f f g , namely 



BL / [f g (u,v) + g(Vu,Vv)]dA. (26) 



fi, 



We express the gradient energy in terms of two Cartesian 
mole fractions u and v, 

g{Vu,Vv) =e„|Vu| 2 + e„|Vi;| 2 + e m ,(Vu- Vu), (27) 

where e u = (3/8)(£ 2 + I 2 ), e v = (l/8)(£ 2 + I 2 ) + (1/2)4 
and e uv — (\/3/4)(^ 2 — For isotropic gradient energy, 



£i = £, we have e u — e v — (3/4)£ 2 and e uv = 0. According 
to the scaling (|23[) from (u,v) to Li, we can express Cl xs 
as 



^d*L 2 I?+g(Vu,Vv) 



i=l 
3 



dA 



dfi^-lf + f |V^| 2 



(28) 



where 

l\ = - 



'e u ui 2 u 31 + e v vi 2 v 31 + ^-(uuvu + V12U31)] 



ft 4 
2 



<-2 - -T4 ^ U U 23 U 12 + £ W U23'fl2 + — ("23^12 + V 23 U 12 )] 
^3 = -^[e«W3l"23 + e v V 31 V 2 3 + ^J-{u 3 lV23 + U31U23)]- 



(29) 



1. Case: di = 1 

Consider a special case of the potential f g in (fTU| when 
all of the weighting coefficients are equal to one, i.e. 
di = 1. According to the asymptotic analysis in our 
previous work [3, sec. II.C], we know how to obtain ana- 
lytical solutions for two-phase transitions in a region far 
from the three-phase contact line for the form (|28|) of the 
excess grand potential f2 xs . These solutions show that 



the interfacial tension <7 Q( 3 is proportional to y i\ + 1\ , 
and the similar relations applied to a ia and cjg 7 (see 0, 
eq. 25]). So we can express the relation ([2]) of the dihedral 
angles and the interfacial tensions as following: 



sin 9p 



suit 



ti + ei 



i\ + i 2 



(30) 



in which the denominators can be obtained by calculating 
the sums of pairs of i 2 in ([29]) . 



i\ + i\ = jj[e u ul 2 + e v v{ 2 + e uv ui 2 vi 2 } 
2 

^2 + ^3 = j^i e uul 3 + e v v 2 3 + e uv u 23 v 23 \ (31) 
2 

^3 + h = T4 [ e « u 31 + e v V 31 + £ U vU 3 iV 31 ] . 



2. Case: di = 1 and ii = . 



Following by the discussion in Sec. IIII C"2l for the spe- 
cial case with equal weighting coefficients di — 1 of the 
potential f g in (|10p . a subset of this special case posses 
an interesting property. This subset is specified by as- 
suming the gradient energy g in (|27p is isotropic, given 
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by £j = I. Then, the sums of pairs of If in (|31l) become 



V 2 



12 



n 



2h 4 
U 2 
' ~ 21? l " 2:; 
U 2 



2h 4 ' 



*31 



y 12j 



3^% 2 
2^ 



(32) 



3^ 
2h 4 



o2 
J 2i 



where Si is the length of the side of the small triangle 
opposite to the vertex (v,i,Vi), defined earlier in the ex- 
pression (fhf]) for the eigenvalues A± of the Hessian matrix 
of f g . Then, the relation of dihedral angles in (|30|) re- 
duces to 



sin 9 a 
Si 



sin f 



sinf 



So 



(33) 



Therefore, for this special case, the sine of each dihedral 
angle of a bulk phase is proportional to the length of the 
side of the small triangle opposite to the vertex corre- 
sponding to the given bulk phase. This is illustrated in 
Fig. [8] and will be revisited in Sec. IIV Al 

Notice that, in general, ^ / ^ ^ ^3 in (|32|) . which 
leads to the lengths Si of the three edges of the small 
triangle unequal. According to the relation (1531) of the 
dihedral angles and Si, the three dihedral angles are dif- 
ferent. Therefore, f g is capable of representing an asym- 
metric three-phase contact line with isotropic gradient 
energy 



Following the variational method introduced in 0, 
sec. II. C] (details in [24], sec. 5.2.3]), we obtain the far- 
field solution for u at the a/3-interface, 



u{s) 



Ul + u 2 , U1—U2 



tanh 



(35) 



where we choose s = as u — (u\ + u 2 )/2 and define the 
interfacial width parameter of the a/3-interface as 



'int, aft 



"12 



(dj + 4) 



V2\ 



d 2 + d 2 2 



(36) 



Alternatively, if u\ 2 = 0, we can use the relation of L3 
to obtain v(s) for the a/3-interface 



vis) = tann 



(37) 



Fig. [7] illustrates the analytic far-field solutions for the 
mole fraction Yi at the a/3-interface for the special case in 
which all of the weighting coefficient are constants d% = 1 
and the gradient energy is isotropic li — 1 (e„ = e v = 3/4 
and e uv = 0). Compared to the similar analytic far-field 
solutions for our original potential f° given by ® in P, 
fig. 5], the mole fraction Y3 is no longer a constant and, 
in general, the solutions for Y\ and Y 2 are not symmetric 
with respect to the interface. 



D. Asymptotic far- field solutions 

Alternatively, we can find the relationship of the di- 
hedral angles of a system specified by the potential f g 
in (|10p directly. We consider a transition from phase 
/9 = {u2,v 2 ) to phase a — (u±,vi), which follows a 
straight line L3 = within the Gibbs triangle as illus- 
trated in Fig. 3(a) In the far-field limit, the transition 
occurs in a regime far from the three-phase contact line 
in the physical space. In this regime, the three mole 
fractions Yi satisfy the boundary condition Vl^ ■ h = 0. 
By assuming ui 2 ^= 0, we can interchange v with u ac- 
cording to L3 = 0, i.e. v — v 2 = (vi 2 /ui 2 )(u — u 2 ) and 
v — Vi — (vi 2 /ui 2 )(u — ui). Then the problem is es- 
sentially a one-dimensional problem. The excess grand 
potential fl xs in (|28f reduces to the form 



h 

U12 



H(u) 



duY 
ds~) 



ds, (34) 



where s is a coordinate perpendicular to the a/3-interface 
measured from (5 to a and w is the width of an area in the 
far field regime as illustrated in the bottom of Fig. [T] (for 
details, see 0, fig. 4]), H(u) = (d\ +d?,)(u-u 1 ) 2 (u-u 2 ) 2 , 
&nde uv = (ul 2 /h 4 )(e u ul 2 +e v v 2 2 + e uv ui 2 vi 2 ). The limits 
of integration are effectively from —00 to 00. 




FIG. 7. Analytic far-field solutions for the mole fraction Yi for 
a special case of the potential f g in (|10|l at the a/3-interface. In 
this special case, di — 1 and li — 1 (e„ = e v =3/4 and e uv = 
0). The three minima of f g are (m,vi) = (0.3087,0.1667), 
(U2,v 2 ) = (0.9060,0.1167), and (u 3 ,v 3 ) = (0.4974,0.6667). 

The interfacial tension of the a/3-interface in the far- 
field limit is given by 



V(dl + d 2 )i uv _ Bh4 y/W + + 



3\u 12 



3a/2 



Bh\ 
(38) 
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We can calculate the interfacial tensions of the ^7- 
intcrfacc and the 7a-interface by the same method. Com- 
pared to the previous result ([30jl of the equilibrium dihe- 
dral angles for a system with potential f g , a more general 
relation obeys 



sin 9 a 



sin Op 



(4 



4m 



(I) 



(d 2 + d 2 )(i 2 + i 2 ) 



sm( 



(39) 



IV. 



GEOMETRY OF INTERFACIAL TENSION 
AND LINE TENSION 



Recall the discussion of the mapping from the physical 
space to the Gibbs space in the Sec. II Bl In a ternary fluid 
system with three phases, the interfacial tensions associ- 
ated with the three interfaces correspond to curved lines 
(trajectories) within the Gibbs triangle, that connect the 
three minima of a given potential. Each trajectory can 
be obtained by minimizing the excess grand potential for 
a two-phase transition in a far-field limit. These trajec- 
tories bound a region of the Gibbs space that associates 
with three-phase transitions. Consequently, the line ten- 
sion associated with the three-phase contact line is deter- 
mined by the excess grand potential within this region. 
This mapping implies a geometrical relationship among 
the values of interfacial and line tensions and the size and 
shape of the core area. 

Here, we demonstrate this geometrical connection by 
some special cases of the specific potential f g given by 
(fTOf , in which the core area is the small triangle formed 
by the minima of the potential within the Gibbs triangle. 



Case: di = 1 and li 



As shown in the Sec. lIII Al for the potential f g in (fTUl) . 
the trajectory for a two-phase transition follows a straight 
line. In the special case of f g discussed in Sec. lIII C~2l all 
of the weighting coefficients di = 1 and the gradient en- 
ergy is isotropic, li = I. Compared to the form ([2J of 
the classical result, in which the sine of a dihedral angle 
is proportional to a corresponding interfacial tension, the 
relation in (|33|) indicates that the sine of a dihedral angle 
is proportional to the length of the side of the small tri- 
angle, within the Gibbs triangle, opposite to the vertex 
for the corresponding bulk phase. This leads to the fact 
that the three interfacial tensions are proportional to the 
three lengths of the sides of the small triangle. It also 
implies that the small triangle is a Neumann triangle, 
which is similar to the computational boundary of the 
physical domain (see Fig. 1(a) and Fig. [8]). 



To prove this property explicitly, we substitute l\ + l\ 



3(y) a-e r 



fS,xfr v 



la 



S.oc a , 



it- 0„ 



FIG. 8. Geometry of the three interfacial tensions in the 
Gibbs space for a special case of the potential f g given by (|10|) . 
in which the weighting coefficients are equal to one di — 1 and 
the gradient energy is isotropic li = 1. The three vertices 1, 2, 
and 3 correspond to the three minima a, ft, and 7 of f g . The 
length of the side of the small triangle opposite to the vertex 
i is labeled by Si, and the area of the small triangle is A. The 
interfacial tension of the a/3-interface a a p is proportional to 
the length of the side of the small triangle that connects a and 
/3, i.e. S3. Similarly, ap-, oc Si and ^70 oc S2. a , 6@, and 9~ f 
are the three dihedral angles corresponding to the three-phase 
contact line. 



from 435J) into ^ to obtain 



JDl 2 

a a R = —^h S3. 
V6 



(40) 



Similarly, for the /37-interface and 7a-interfaces, the ten- 
sions are 



Bl. 



a B~t = — -^h^Si and <7~ Q = — r=h 2 S2- 
V6 7 V6 



B. Case: d x = 1, & = £, and S t = S 



(41) 



We consider a subset of the case treated in Sec. IIV A[ 
in which the small triangle formed by the three minima 
of the potential f g in (fl0|) is equilateral, i.e. Si = S. 
Then, the area of the small triangle is 



2 4 



(42) 



From the relation (|3"2"|) between the sums of pairs of £| 
and Si, we find that 



|a _ p _ 3£ 2 S 2 U 2 h yfll 2 



2h 4 2 2h 4 y/3 2h 3 



(43) 



According to the form (l2"rJl) of the excess grand poten- 
tial fl xs , we use the Kerins-Boiteux formula [26j to obtain 
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a line tension associated with the three-phase contact 
line. 



t = B [-f g (u,v)+g{Vu,Vv)]dA 

J A 

3 



-ZfiZi-lf + jlVZi] 2 



dA 



(44) 



where f is the dimensionless quantity introduced in [1], 
eq. 53]. Then, from the expressions for area A in 
and l 2 in 03]), we find 



(45) 



which is proportional to the area of the small triangle, as 
illustrated in Fig. [9] 




FIG. 9. Geometry of the line tension in the Gibbs triangle for 
a special case of f g with equal weighting coefficients di = 1, 
isotropic gradient energy £i — £, and the small triangle formed 
by the three minima is equilateral. The line tension and the 
interfacial tensions are proportional to the area A and the 
lengths of the sides of the small triangle, respectively. 

From above cases, the small triangle in the Gibbs space 
not only tells us the compositions of the three bulk phases 
in terms of mole fractions from the locations of its three 
vertices; its size and shape also give us information about 
the equilibrium dihedral angles and the relative strengths 
of the interfacial tensions. Moreover, the line tension is 
proportional to the area of this small triangle. 



sum of three quartic polynomials, each of which vanishes 
along two parallel lines in the Gibbs triangle. The three 
minima of our potential are located at the intersections of 
any three among these six lines. By recognizing this geo- 
metrical structure, we can construct generalized quartic 
potentials with three minima arbitrarily located at the 
vertices of a small triangle within the Gibbs triangle. 

As a first generalization, for each vertex, we can choose 
a linear function that vanishes along a line that passes 
through a pair of vertices and another linear function that 
vanishes along a parallel line that passes through the re- 
maining vertex. The squares of these two linear functions 
form paraboloids that vanish at these two lines. We use 
the product of these two paraboloids to construct one 
of the three quartic polynomials whose sum is the gen- 
eralized potential. By performing similar constructions 
at the other two vertex pairs, we can build a potential 
which is the sum of three quartic polynomials that van- 
ishes at only three points. Our original potential belongs 
to a subset of this generalization where the two "vanish- 
ing lines" of each quartic polynomials are parallel. To 
construct a potential in this subset, we need six param- 
eters, namely the coordinates of the three minima and 
three weighting coefficients for each quartic polynomi- 
als. However, we find that each paraboloid that vanishes 
at a line passing through one vertex can be generalized 
to another paraboloid that only vanishes at a point, the 
vertex. By means of this second generalization, we con- 
struct another potential with nine parameters related to 
the shapes of the new paraboloids, and containing six 
parameters to characterize the three minima. We check 
the properties of these generalized potentials by studying 
their first and second derivatives. 

Generalized potentials in which the valleys connecting 
any of the two minima follow straight lines in the Gibbs 
triangle have simple analytic far-field solutions. We prove 
that this subset of generalized potentials is a special case 
of our first generalization that was constructed by three 
pairs of parallel lines. Analytic far-field solutions of this 
subset are presented. Moreover, we find that this subset 
is a linear transformation of our original potential. By 
scaling, we can relate their solutions to those for our orig- 
inal potential. When the weighting coefficients of the po- 
tentials in this subset are equal and the gradient energy 
is isotropic, the lengths of the sides of the small trian- 
gle formed by the three minima are proportional to the 
corresponding interfacial tensions. For the case of equal 
interfacial tensions, we are able to calculate a line tension 
that is proportional to the area of the small triangle. 



SUMMARY AND CONCLUSIONS 
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Appendix A 

We can generalize the potential f g in (flTJ)) by replacing 
Ii with an arbitrary linear function Ji that passes through 
(Ui, Vi). If Ji(u, v) = CLi{u — Uj) + bi{v — Vi), where a* and 
hi are arbitrary constants, we can obtain a generalized 
potential 



h = 22Li(u,v) 2 Ji(u,v) 2 , 



i=l 



(Al) 



and 



as illustrated in Fig. 10(a) In this case, Ji 
Li — do not need to be parallel. 

For the special case in which J\ , J 2 , and J 3 match with 
L 2 , L 3 , and Li, we can define a potential that employs 
only three lines, namely 



f 2 =d\ 2 Li(u,v) 2 L 2 {u,v) 2 
+ dl 3 L 2 (u,v) 2 L 3 (u,v) 2 



dliL 3 (u,v) 2 ii(u,w) 2 , 
(A2) 



where di 2j d 23 , and d 3 \ are weighting coefficients. At each 



minimum, only two lines meet as illustrated in Fig. 11(a) 
The contours of potential f 2 are shown in Fig. ll(b)| 



Note that /1 in (|A1[) and f 2 in (|A2|) are special cases 
of the potential fo defined in (fT7|) . 



Appendix B 

To study the general form of the subset of general- 
ized potentials in which any two-phase transition follows 
a straight line, we explore the transition between two 
phases by means of a potential f 2 n that is a quartic func- 
tion of three mole fractions having three minima (u±, V\), 
(u 2 ,v 2 ), and (u 3 ,v 3 ) located at the three vertices of the 
Gibbs triangle. The potential function along any two 
of the three wells actually reduces to a one-dimensional 
potential /id- Consider the transition from (ui,vi) to 
(u 2l v 2 ) that follows a straight line L 3 — 0. We let 



/2£>li„ =0 ( u > v ) = f- 



1D- 



(Bl) 



where /id is a positive quartic potential function having 
two minima at (ui,vi) and {u 2 ,v 2 ). Since fm follows a 
straight line L 3 = within the Gibbs triangle, we can 
express /id as a function of one variable. Then, we let 



Iid{u) = (u- ui)(u - u 2 )g(u), 



(B2) 



where g(u) is a quadratic function of u. By requiring 
dfm/du = at ui and u 2 and 9 2 /id /du 2 > 0, we obtain 



3(y), 



,L, = 



la 



/, = 



J,=0. 



{p) 



(a) Geometry of the three minima of /1 




(b) Contours of /1 

FIG. 10. |(a)| Geometry of the three minima of the quartic 
potential Ji given by (|A1|) . a — (ui,Vi), f3 = (u2,«2), and 
7 = (113,1)3) are the three minima of this potential arbitrar- 
ily located at any three points in Gibbs space. Li — rep- 
resents a line that passes through (v,j,Vj) and (uk,Vk) for 
j 7^ k 7^ i. (b) Contour plot of potential /1 for (ui,ui) 



= (0.3087,0.1667), {ui,vi) = (0.3087,0.1667), (u 2 ,v 2 ) = 
(0.9060,0.1167), (u s ,v 3 ) = (0.4974,0.6667), and (01,61) = 
(1,0.1), (o 2 ,6 2 ) = (-1,0.1), and (03,63) = (-0.1,1). 



We can generalize the contribution of /id(u) to 

f 2 n(u,v), which is denoted as f^r> by combination of 
Li(u,v) and K\{u, v), i.e. 



J 2D 



.^(•u, u)isTi(ii,u), 



(B4) 



similar to the potential fa in pT|) , which vanishes at the 
three minima. However, if the transition between (u\, v\) 
and (^2,^2) follows the straight line L 3 — 0, we need 



dLf(u, v)Ki(u, v) 



(B5) 



L 3 =0 



/id(m) oc (u - ui) 2 (u - u 2 y 



(B3) where L 3 = {v\ 2 , —u\ 2 ) is a normal vector perpendicular 
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3(y), 



^ L,=0 



1(a) 



L, = 



2(P) 



to L3 = in the Gibbs space. Thus, we calculate 



d{L\K x ) 



dL 3 



\L 3 =0 



= 2 



V12 



{v - vi)(v - v 2 ) 



x { {v 12 f 



ai h 01 + Ci 

«12 



-"12 



"12 



7 "12 1 2 2 

ai&i h 0, + c, 

U12 



(v - v 2 ) 

(B6) 



(a) Geometry of the three minima of f 2 




(b) Contours of fa 



FIG. 11. (a) Geometry of the three minima of the quartic 
potential fa given by (|A2I) . a = /3 — (112,^2), and 

7 = (113, v$) are the three minima of this potential arbitrar- 
ily located at any three points in Gibbs space. Li = rep- 
resents a line that passes through (uj,Vj) and (ufe,Wfe) for 
j 7^ k 7^ i. (b) Contour plot of potential /2 for (ui,Vi) 
= (0.3087,0.1667), (wi,«i) = (0.3087,0.1667), (1*2,-02) = 
(0.9060,0.1167), («3,«3) = (0.4974,0.6667), and d 12 = d 23 = 
d 3 i = 1. 



in which we substitute (u — u\) — {uw/vvijiy — i>i) for 
L3 = by assuming vn 7^ 0. In general, h ^ 0, U12 7^ 0, 
(« — «i) 7^ 0, and (u — U2) 7^ 0. The nontrivial solution 



for (BSD is «i 
reduces to 



-U12, and ci = 0. Then, 



(B7) 



Similarly, we can apply this same analysis to the transi- 
tions of (u2, V2) (1*3, U3) and (1*3, 1)3) f-> (ui, «i). Thus, 
we conclude that 



flD — fg- 



(B8) 



In Section |TT1 we have shown that the potential f g has 
our desired properties. Moreover, the two-phase transi- 
tions of f g follow straight lines in the space of the Gibbs 
triangle, which allows us to find analytical solutions in 
the far-field. 



* changyoul@gmail.com; Corresponding author 

* widom@andrew.cmu.edu 

* sekerka@cmu.edu 

[1] C.-Y. Lin, M. Widom, and R. F. Sekerka, Physical Re- 
view E 85, 011120 (2012) 

[2] J. W. Gibbs, The Collected Works of J. Willard Gibbs: 
Thermodynamics, Vol. 1 (Longmans, Green, 1928) 

[3] J. S. Rowlinson and B. Widom, Molecular theory of cap- 
illarity (Dover, 2002) 

[4] F. E. Neumann, "Vorlesungen fiber die theorie der cap- 
illaritat (ed. a. wangerin)," (Teubner, Leipzig, 1894) 
Chap. 6 §1, pp. 161-2 

[5] A. W. Adamson and A. P. Gast, Physical chemistry of 
surfaces (Wiley-Interscience, 1997) 



[6] J. Gaydos and A. W. Neumann, Journal of colloid and 

interface science 120, 76 (1987) 
[7] J. Drelich, J. D. Miller, and J. Hupka, Journal of colloid 

and interface science 155, 379 (1993) 
[8] M. Grunze, Science 283, 41 (1999) 
[9] B. H. Weigl and P. Yager, Science 283, 346 (1999) 
[10] A. I. Hienola, P. M. Winkler, P. E. Wagner, 

H. Vehkamaki, A. Lauri, I. Napari, and M. Kulmala, The 

Journal of chemical physics 126, 094705 (2007) 
[11] E. Sackmann and R. F. Bruinsma, ChemPhysChem 3, 

262 (2002) 

[12] J. Fukai, Y. Shiiba, T. Yamamoto, O. Miyatake, 
D. Poulikakos, C. M. Megaridis, and Z. Zhao, Physics 
of Fluids 7, 236 (1995) 



14 



[13] D. Bonn, J. Eggers, J. Indekeu, J. Meunier, and E. Rolley, 
Reviews of modern physics 81, 739 (2009) 

[14] J. O. Indekeu, Physica A: Statistical Mechanics and its 
Applications 389, 4332 (2010) 

[15] A. Amirfazli and A. W. Neumann, Advances in colloid 
and interface science 110, 121 (2004) 

[16] L. Schimmele, M. Napiorkowski, and S. Dietrich, The 
Journal of chemical physics 127, 164715 (2007) 

[17] D. M. Anderson, G. B. McFadden, and A. A. Wheeler, 
Annual review of fluid mechanics 30, 139 (1998) 

[18] L. D. Landau and E. M. Lifshitz, Phys. Zeit. Sowjetu- 
nion 8 (1935), see also, "Electrodynamics of Continu- 
ous Media", §45, p. 158, reprinted by Beijing World Pub- 
lishing Cooperation by arrangement with Butterworth- 
Heinemann (1999) 

[19] A. Umantsev, Field Theoretic Method in Phase Transfor- 



mations, Vol. 840 (Springer New York, 2012) 
[20] I. Szleifer and B. Widom, Mol. Phys. 75, 925 (1992) 
[21] K. Koga, J. O. Indekeu, and B. Widom, Faraday Discuss. 

146, 217 (2010) 
[22] K. Koga and B. Widom, J. Chem. Phys. 128, 114716 

(2008) 

[23] K. Koga, J. O. Indekeu, and B. Widom, Physical review 
letters 104, 36101 (2010) 

[24] C. Lin, Mean-field density functional theory of a three- 
phase contact line (PhD thesis, Carnegie Mellon Univer- 
sity, 2012) 

[25] C. Varea and A. Robledo, Physica A: Statistical and The- 
oretical Physics 183, 12 (1992) 
[26] J. Kerins and M. Boiteux, Physica A 117, 575 (1983) 



